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Abstract. An improved cloud/snow screening technique in 
the Multi- Angle Implementation of Atmospheric Correction 
(MAIAC) algorithm is described. It is implemented as part of 
MAIAC aerosol retrievals based on analysis of spectral resid- 
uals and spatial variability. Comparisons with AERONET 
aerosol observations and a large-scale MODIS data analy- 
sis show strong suppression of aerosol optical thickness out- 
liers due to unresolved clouds and snow. At the same time, 
the developed filter does not reduce the aerosol retrieval ca- 
pability at high 1 km resolution in strongly inhomogeneous 
environments, such as near centers of the active fires. De- 
spite significant improvement, the optical depth outliers in 
high spatial resolution data are and will remain the prob- 
lem to be addressed by the application-dependent specialized 
filtering techniques. 


1 Introduction 

A new Multi- Angle Implementation of Atmospheric Correc- 
tion (MAIAC) algorithm developed for Moderate Resolution 

Imaging Spectroradiometer (MODIS) was described recently 

(Lyapustin et al., 201 la, b). This is a generic algorithm which 
retrieves aerosol information over land simultaneously with 
parameters of the bidirectional reflectance distribution func- 
tion (BRDF) model. To achieve this goal, MAIAC uses the 
time series of MODIS radiance measurements as well as pro- 
cessing of groups of pixels. This approach utilizes the dif- 
ference in the time- space variability of aerosols and surface 
reflectance which can be captured with the daily near-global 


coverage of MODIS: namely, aerosols vary slowly in space 
but may change between consecutive MODIS observations, 
whereas the land surface reflectance has a high spatial vari- 
ability but low rate of change at short time intervals. A sim- 
ilar idea has recently been implemented for the advanced 
processing of PARASOL data (Dubovik et al., 2011). 

MAIAC aerosol retrievals are performed at high 1 km res- 
olution which is needed in different applications such as vis- 
ibility assessments (Wang et al., 2009), aerosol source iden- 
tification, air quality analysis (Hoff and Christopher, 2009) 
etc. In a recent work, Emili et al. (2011) evaluated MAIAC 
cloud/snow mask and aerosol products in the region of Euro- 
pean Alps characterized by a heterogeneous aerosol distribu- 
tion with strong impact of topography and aerosol sources 
localized in the narrow valleys with width of several km. 
While this study clearly demonstrated benefits of the high 
resolution data as compared to the standard 10 km MODIS 
product (Levy et al., 2007), including improved spatial cov- 
erage and 50 % increase in the number of observations, it 
has also revealed residual cloud and snow contamination. 
This effect becomes particularly noticeable in rather pris- 
tine Alpine conditions with low average mid- visible aerosol 
optical thickness AOT ~ 0.05-0.2. The problem of bias was 
successfully overcome by Emili et al. (2011) with AOT data 
filtering where the main filter was based on the 3 x 3 pixel 
spatial variance test (a < 0.05). In more detail, this filter suc- 
cessively removed the highest AOT value from the 3x3 km 2 
area if the standard deviation exceeded 0.05, and then aver- 
aged the remaining values effectively leading to 3 km res- 
olution of the aerosol product. The filtering significantly 
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improved correlation of MAIAC data with AErosol RObotic 
NETwork (AERONET) AOT for the selected mountainous 
sites, for example from 0.2 to ~0.8 for Laegeren and 
Davos, Switzerland, by excluding ~30 % and ~50 % of AOT 
retrievals, respectively. 

While filtering is clearly required for some applications, 
such as climatology analysis, it would have negative con- 
sequences for the others. For example, the a -filter along 
with reduction of the effective spatial resolution of MAIAC 
AOT from 1 km to 3 km would eliminate many meaning- 
ful retrievals with point sources such as fire smoke plumes. 
It is, therefore, desirable to address the problem of resid- 
ual cloud/snow contamination within MAIAC itself. In this 
work, we explore the opportunities which exist within MA- 
IAC aerosol algorithm to improve the cloud/snow mask. 


2 Spectral residuals and spatial variability 

The full description of MAIAC has been given before (Lya- 
pustin and Wang, 2009; Lyapustin et al., 2011a, b). Below, 
we provide the minimum level of detail which is only rel- 
evant for the current discussion. MAIAC processing starts 
with gridding MODIS LIB data to a 1 km resolution (Wolfe 
et al., 1998). The gridded data are placed in the Queue 
which stores from 5 (poles) to 16 (equator) days of im- 
agery, depending on latitude. The Queue implements a slid- 
ing window algorithm used for cloud masking (CM) and 
surface characterization. Both algorithms utilize 1 km grid 
cells, which are called pixels, as well as fixed 25 x 25 km 2 
areas called blocks. MAIAC CM algorithm (Lyapustin et al., 
2008) provides a generally robust performance which is sim- 
ilar to that of the MODIS operational cloud mask (Ackerman 
et al., 1998) or may exceed it in difficult conditions, for ex- 
ample over bright surfaces and snow. Common to all CM 
algorithms, the MAIAC CM has a limited ability to identify 
thin or sub-pixel clouds. MAIAC CM algorithm includes a 
dynamic land- water- snow classification based on the time se- 
ries analysis. The snow detection tests are commonly based 
on the fixed thresholds, which automatically creates a prob- 
lem of residual snow contamination in aerosol retrievals. The 
main outcome of the MAIAC surface characterization, rel- 
evant for aerosol retrievals, are the BRDF model parame- 
ters and surface reflectance uncertainty (£x) at the top of 
atmosphere (TOA) for every 1 km grid cell in the reflective 
MODIS bands. 

The aerosols are modeled conventionally as a superposi- 
tion of the fine and coarse modes. Following the MODIS op- 
erational Dark Target algorithm MOD04 (Levy et al., 2007), 
the fine and coarse aerosol models in MAIAC are fixed re- 
gionally based on AERONET (Holben et al., 1998) climatol- 
ogy. MAIAC uses the latest MODIS measurements to per- 
form aerosol retrieval based on the knowledge of spectral 
surface BRDF and its uncertainty in bands B3 (0.47 pm), B1 
(0.67 pm) and B7 (2.13 pm) from the previous retrievals. For 


each pixel (/, j), it computes AOT by matching the mod- 
eled TOA reflectance to the measurement in the Blue band 
(B3). This procedure is repeated for the increasing values of 
the coarse mode fraction (CMF) characterized by parameter 
rj equal to ratio of volumetric concentrations of the coarse 
and fine fractions. The final solution (r, ij) is selected based 
on the rmse test which is computed using the Red (Bl) and 
shortwave infrared (SWIR, B7) bands: 

1 ( ^ eas ' X - Rjp°' (r k (r])) 1 2 

— i — I 

< 1 or x,j(r]) = min. (1) 

If Eq. (1) cannot be satisfied with aerosol models, the algo- 
rithm also tries a liquid water cloud model. The latter rep- 
resents a cloud consisting of 5 pm water droplets with nar- 
row size distribution (o = 0.5 pm) which is used to test pos- 
sible cloud contamination in each pixel. This additional test 
improves cloud detection capturing many thin or small sub- 
pixel clouds (e.g., see Fig. 5 from Lyapustin et al., 2011b). 
Prior to aerosol retrievals, a snow test (Li et al., 2005) origi- 
nally implemented in the MOD04 algorithm is performed to 
filter undetected snow pixels. 

While the rmse (x) test proved to be useful for improved 
cloud masking, there is additional information contained in 

( y^Meas _ y^Thcor \ 

— 1. For 

example, a retrieval for a thin cloud pixel with the back- 
ground aerosol model will result in positive residuals <$o. 67 , 
<$ 2 . 13 . For a partly cloudy pixel, the residuals will be positive 
with the aerosol models and will change sign when the cloud 
model is used in the retrievals. 

2.1 Spectral residuals 

A proposed simple cloud test is based on the difference in 
spectral dependence of extinction of aerosols and clouds due 
to a large difference in the particle size. To understand its ca- 
pabilities and assess sensitivity to thin clouds over different 
surfaces, numerical simulations were conducted. The TOA 
radiance was first simulated for a given atmosphere- surface 
model using code SHARM (Lyapustin, 2005), and then MA- 
IAC aerosol retrieval was applied. We used two surface types 
representing a typical dense vegetation and a bright urban 
area whose BRDF model and its uncertainty were provided 
by MAIAC from MODIS data. The green and bright ur- 
ban areas geographically represent the summertime north- 
ern Washington DC, with albedo g = {0.014, 0.02, 0.033, 
0.061} and g = {0.04, 0.081, 0.149, 0.20}, respectively, as a 
measure of surface brightness in the MODIS channels B8 
(0.412pm), B3, Bl and B7. Here, the MODIS “Deep Blue” 
band B8 was added to the set of channels used by MAIAC 
in aerosol retrievals. The liquid water cloud was modeled us- 
ing a lognormal size distribution with radius 10 pm and stan- 
dard deviation 0.5 pm using refractive indices from Hale and 
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Fig. 1. Spectral residuals (60 from aerosol retrieval for a cloudy pixel with optical thickness of 0.234 (a-b) and 0.7 (c) at MODIS wavelengths 
0.412, 0.47, 0.67 and 2.13 pm. The results are shown for two surface types representing dense vegetation ((a), GP) and an urban (city) area 
((b-c), BP). The red, black and blue lines correspond to cosines of view zenith angle fi = — 1, —0.7 and —0.4, and the solid, dashed and 
dotted lines represent relative azimuths of 35°, 90°, and 145°, respectively. 


Querry (1973). A dynamic East Coast aerosol model (see 
Lyapustin et al., 201 lb) was used in the retrievals. 

The test results are presented in Fig. 1. The plots (a-b) 
show spectral residuals for the vegetated and urban surfaces 
obtained with aerosol model for a thin cloud pixel with op- 
tical thickness (COT) of 0.234, and plot (c) shows results 
for the urban surface and a thicker cloud of COT = 0.7. The 
different lines represent different view geometries: the red, 
black and blue lines correspond to cosines of view zenith 
angle /x = cos(VZA) = — 1 , —0.7 and —0.4, while the solid, 
dashed and dotted lines represent three relative azimuths of 
35° (forward scattering), 90°, and 145° (backscattering). The 
residual in the Blue channel (0.47 pm), which is used to 
compute AOT, is always zero. It is positive in the Red and 
SWIR bands over the dark and vegetated surfaces. The plot 
(a) shows that a simple criterion 6o.67 > 1.5-2, 62.13 > 1.5- 
2 could be used in this case to detect very thin liquid wa- 
ter clouds with COT ~ 0.25. Over brighter surfaces, however, 
the residuals may take both positive and negative values, de- 
pending on the view geometry. In this case, a sufficient de- 
tection sensitivity is attained for thicker clouds (COT > 0.7), 
as illustrated by plot (c). These plots also show that adding 
the “Deep Blue” channel (0.412 pm), where the residual sys- 
tematically takes negative values, may enhance the cloud dis- 
crimination capability of the proposed spectral test. 

A similar idea can be used for detection of the residual 
snow which increases surface brightness in the visible wave- 
lengths (60. 67 > 0) and decreases it in the SWIR (62.13 < 0). 

While the idea of the proposed spectral test is seem- 
ingly simple, its realization is complicated by several fac- 
tors: (1) the spectral surface reflectance in MAIAC is known 
with uncertainty characterized by its standard deviation e^. 
This error includes contributions from all sources including 
gridding, atmospheric correction and fitting limitations of the 
BRDF model. Assuming Gaussian distribution of errors, the 
specific reflectance would agree with the model to within 
±sx in ~68 % cases and within ±2 ex in ~95 % cases; (2) the 
surface can change since the last BRDF retrieval. For exam- 
ple, rain can darken the soil decreasing its reflectance in both 
Red and SWIR channels, whereas undetected sub-pixel snow 


would increase surface reflectance in the Red band and de- 
crease it in the SWIR. Over vegetated surfaces, the common 
perennial changes are related to the vegetation phenology 
tracking transitions between winter and summer in the north- 
ern latitudes or between wet and dry seasons in sub-tropics. 
While the “green-up” surface signal is spectrally unique, the 
senescence or “browning” of the surface presents a particu- 
lar problem because it is spectrally similar to the effect of 
thin clouds in the Red-SWIR bands. This discussion shows 
that the individual pixel tests are prone to errors, and they 
should be used together with the larger- scale analysis based 
on groups of pixels (blocks), which would provide statistical 
mitigation of errors and a robust separation between clouds 
and surface change. 

2.2 Aerosol spatial variability 

With the block-level analysis, justified above, one can use 
additional spatial variability techniques to screen outliers 
caused by clouds and snow. They are based on a relative ho- 
mogeneity of global aerosols at scales below ~50km (e.g., 
Anderson et al., 2003). For example, the MOD04 Collec- 
tion 5 algorithm uses a 3 x 3 spatial variance cloud filter 
(Martins et al., 2002) and discards the darkest 20% and 
brightest 50% pixels in the 10 x 10 km 2 box, helping to 
screen cloud shadows and clouds/snow, respectively. A sim- 
ilar approach was implemented in MAIAC validation anal- 
ysis against AERONET (Lyapustin et al., 2011b, c) based 
on screening of the high 50% of retrieved AOT data. For 
this reason, MAIAC validation was not generally affected by 
the outliers. In addition, averaging the remaining data over a 
10-20 km window allowed us to account for the meteorolog- 
ical conditions and the time difference between AERONET 
measurements and MODIS overpass (Ichoku et al., 2002), 
as well as to increase the comparison statistics. Emili et 
al. (2011), however, conducted validation using a single 1 km 
AOT value closest to the AERONET sunphotometer location 
which explicitly revealed the high scatter and biases in the 
unfiltered MAIAC AOT data over mountainous regions. 
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Fig. 2. Block-diagram of MAIAC aerosol algorithm with enhanced 
cloud- snow- shadow masking. The symbol & stands for ‘AND”. 


The spatial variability analysis, based on the 25 km blocks, 
was introduced in the current version of MAIAC. Specifi- 
cally, it filters high AOT values when clouds and/or snow 
are detected by the CM algorithm. The threshold linearly de- 
pends on the cloud fraction (CF), decreasing from 60th per- 
centile for CF = 0.05 down to 20th percentile for CF = 0.7. An 
extensive analysis of MODIS data showed that the dynamic 
threshold depending on cloud fraction provides much better 
cloud screening than the static global threshold. If the cloud 
fraction exceeds 70 %, MAIAC does not perform processing 
for the given block. In cloud-free conditions with snow de- 
tected, the high threshold represents the 25th percentile of 
AOT data. 

The described screening is not applied when the cloud 
fraction is low to preserve MAIAC capability for high 
resolution aerosol retrieval near the aerosol sources. 


3 Algorithm implementation 

The improved scheme of MAIAC aerosol retrievals is illus- 
trated in Fig. 2. It includes the surface change detection com- 
ponent and enhanced cloud, cloud shadow, and snow mask- 
ing. The numbered rectangles on the left represent separate 
routines with the block-level (25 x 25 km 2 ) scope of appli- 
cation. Symbol B on the left indicates that an operation is 
applied to every cloud-free pixel of the block, whereas RB 
means “the Remaining pixels of the Block” unaffected by 
the previous actions. Because the individual pixel tests do not 
guarantee the correct result, the algorithm is designed to al- 


low commission errors, for example detecting clouds in clear 
conditions, with restoration of the correct results in the final 
stage of processing. 

The aerosol algorithm starts with computing aerosol op- 
tical thickness and rmse (r 1 , X 1 )?'; using the background 
aerosol model (k = 1) for every clear pixel of the block ac- 
cording to the MAIAC cloud mask (step 1). A copy of these 
results is saved for a later “Clear Sky Restore” analysis 
(step 9). 

The next retrievals should be repeated with higher CMF 
values (j]) searching for a pair (t, r]) that minimizes the rmse 
given by Eq. (1). At this stage, the undetected surface change 
may introduce a systematic error as was mentioned above. 
For example, during senescence, when the surface is bright- 
ening in the Red and SWIR channels, the straightforward ap- 
proach would overestimate both CMF and AOT and would 
also result in a high commission error of false cloud detec- 
tion. To avoid that, a Surface Change Detection algorithm is 
implemented in step 2. It is applied when AOT is low and 
the day is clear which for a given block is verified by a high 
covariance (cov > HIGH) between the measured reflectance 
(in Bl) and the reference clear-sky image maintained by the 
CM algorithm (Lyapustin et al., 2008). 

The Surface Change Detection looks for a simultaneous 
anti-correlated change of surface reflectance in the Red and 
NIR (band B2) bands during the 16-day interval. To enable 
reflectance comparisons measured at different view angles 
on different days, the surface reflectance is first normalized 
to the standard view geometry of nadir view and solar zenith 
angle 45° using the BRDF model. While the green-up change 
is spectrally unique and can be accepted for an individual 
pixel, the “browning” can be easily confused with undetected 
thin cloud. For this reason, detected “browning” is confirmed 
only if it is observed for at least a quarter of the block’s pix- 
els. Otherwise, the detected “browning” is classified as a ran- 
dom noise and is discarded. The details of this algorithm will 
be described separately. 

If the retrieved AOT is low or rmse < 1 or the surface 
change has been detected, the algorithm reports AOT for 
the background aerosol model in step 3, and aerosol pro- 
cessing for a given pixel terminates. Otherwise, cloud test 
1 (CT1) is applied (<$o.67 > 1.5, 62.13 > 1.5, and 60.412 < 0), 
and if successful, the pixel is flagged as possibly cloudy 
(CM_PCLOUD) . Note that all newly detected CM_PCLOUD 
pixels must be validated in the final “Clear Sky Restore” test 
in step 9. 

Step 4 shows the standard aerosol retrieval loop with 
higher coarse mode fractions (index k) according to Eq. (1). 

The further processing is designed to detect additional 
clouds, cloud shadows and snow. Step 5 helps avoid 
unnecessary processing in clear conditions (cov > HIGH). 

In the next Step 6, the aerosol retrievals are repeated with 
the cloud model for the remaining pixels of the block. The 
pixel is masked as possibly cloudy if rmse k < 1, or if spec- 
tral residuals are negative with the cloud model but were 
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Fig. 3. MAIAC-AERONET AOT 0.47 scatterplots for Laegeren, Switzerland, and Ispra, Italy, based on MODIS Terra 2000-2008 data. The 
four plots in a row represent MAIAC data at original 1 km resolution, and for a 3, 5 and 10 km window. 


positive with the last aerosol model (X(f 67 < 0 , xfn < 0 an d 
XoV > 1, X 2 Y 5 > 1) which often indicates presence of the 
sub-pixel clouds. 

Step 7 implements the spatial variability analysis dis- 
cussed in Sect. 2.2 based on low (AOT 0 . 1 ) and high 
(AOTffigh) AOT thresholds, where the latter depends on the 
cloud fraction. 

The next step 8 performs additional tests for residual snow, 
shadows and clouds: 

- If snow has been detected in the block (N snow >5), then 
the pixel is flagged as CM_CLEAR_SNOW if <5o.67 > 2, 
$2.13 < —2, t ij > AOTnigh and cov is high. If covariance 
is low, indicating presence of clouds, then the pixel is 
flagged as possibly cloudy (CM_PCLOUD). 

- The shadow test <5o.67 < —2, < 52.13 < —2, r ij < AOT 0.1 
provides enhancement to cloud shadows detected by 
MAIAC CM algorithm. 

- An additional cloud test screens pixels with simulta- 
neously high AOT and rmse , Xij > 2, r ij > AOTnigh 
(CT4). 

Earlier we mentioned that uncertainty in the knowledge of 
surface reflectance and surface change often result in selec- 
tion of unrealistically high aerosol coarse mode fraction (and 
high AOT value) or false cloud detection. The last Clear Sky 
Restore Test 9 is designed to correct these errors and restore 
the value of cloud mask and AOT/CMF. It is based on the 
idea that aerosol variability at scale of 25 km is expected to 
be low in clear conditions, so the pixel AOT should be close 
to the average value. To this end, the block-average value 
(AOT av ) is first computed for the pixels with the background 
aerosol model retrieval. Next, for the other cloud-free pixels 


we check if the original AOT value retrieved with the back- 
ground model and saved at stage 1 (r- 1 .) is close to the AOT av . 
Specifically, we restore the CM_CLEAR value of the cloud 
mask for pixels masked as CM_PCLOUD, or the background 
model for pixels with high CMF if ij- < 1.2 x AOT av . 

4 Performance analysis 

The performance of the upgraded algorithm was first as- 
sessed in comparison with AERONET data for different 
sites, including Laegeren, Switzerland, and Ispra, Italy an- 
alyzed by Emili et al. (2011). We used three different averag- 
ing window sizes, 3, 5 and 10 km, and applied a standard val- 
idation approach which filters high MAIAC AOT data above 
70th percentile (Lyapustin et al., 2011b) and requires at least 
3 valid retrievals in the window. In order to test MAIAC AOT 
at 1 km resolution, we used the nearest valid pixel within 
1 km from the AERONET sunphotometer location. 

The results for the Laegeren and Ispra sites, based on 
MODIS Terra 2000-2008 data, are presented in Fig. 3. The 
original 1 km data show the most improvement in correla- 
tion coefficient, from R ~ 0.2 to R ~ 0.84 for Laegeren and 
from R~0.67 to R~0.91 for Ispra, despite the relatively 
high scatter and offset. The number of comparison points is 
the lowest at 1 km increasing by a factor 4-8 for 3 and 10 km 
windows for Laegeren, and 6-14 for Ispra, respectively. The 
regression parameters improve noticeably when switching to 
a 3 km window due to filtering residual outliers and data av- 
eraging, however further improvement with increasing the 
window size to 10 km is only incremental. Moreover, the re- 
gression coefficient for Ispra slightly drops when the window 
size becomes larger than 3 km. Analysis of MAIAC retrievals 
shows a significant heterogeneity of aerosols on hazy days 
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Fig. 4. Illustration of aerosol variability in 50 km area near Ispra 
during 13-17 October 2001. The five columns show the MODIS 
Terra TOA data and MAIAC products including cloud mask, RGB 
NBRF (bidirectional reflectance computed from the BRDF model 
for a fixed view geometry of nadir view and 45° solar zenith angle), 
BRF (or surface reflectance) and AOT 0 . 47 . The blue and light blue 
colors of the cloud mask indicate clear over the land and water, 
respectively, and red-yellow show detected clouds. 

in the area of Ispra created by aerosol transport from highly 
industrialized and polluted Po valley to the south-east and 
blocking effect of Alps on the north-west. Figure 4 gives an 
illustration of significant south-east to north-west AOT gra- 
dient across the 50 km area centered at Ispra during 5 days 
in October 2001. In these circumstances, the smaller 3 km 
window may give a better representation of the local condi- 
tions sampled by the AERONET sunphotometer. A similar 
observation was reported by Lyapustin et al. (201 lc) for Bei- 
jing where a better agreement with AERONET was observed 
with 10 km rather than 20 km window. 

One of the key requirements for the developed cloud/snow 
filter in MAIAC was to keep the aerosol retrievals with high 
spatial variability unaffected, in particular near sources. Fig- 
ure 5a-b gives a large-scale illustration of the algorithm per- 
formance for two days (DOY 235 and 274) of 2003 during 
wild fires in Georgia, USA (e.g., see Christopher et al., 2009 
for a similar event during 2007). Both figures show gridded 
MODIS Terra RGB TOA image, MAIAC RGB BRF, cloud 
mask and aerosol optical thickness for 600 km tiles. The first 
day is characterized by complex meteorological conditions 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 


Fig. 5. (a) Illustration of a large-scale MAIAC algorithm perfor- 
mance during wild fires of 2003 in Georgia (USA) on day 235 (a) 
and 274 (b). The images show MODIS Terra TOA data, MAIAC 
RGB BRF, cloud mask and AOT 0.47 for 600 km tiles. The cloud 
mask color coding is described in caption of Fig. 4. (b) The same as 
Fig. 5a but for day 274 of 2003. 
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with variable cloudiness and high aerosol levels in the middle 
of the image transported by easterly winds. The aerosol con- 
centration is much lower on the northern and southern parts 
of the image where the good quality of cloud masking be- 
comes apparent. The second image (day 274) shows a num- 
ber of active fires and smoke plumes. Despite the high spatial 
aerosol variability, these retrievals generally are not affected 
by the developed filter even in the centers of the burning area. 

5 Conclusions 

Contamination of aerosol retrievals by unresolved clouds and 
snow is a serious problem of remote sensing gaining a par- 
ticular significance at high spatial resolution. Due to a con- 
tinuous spectrum of cloud/snow presence in measurements 
in terms of sub-pixel area and reflectance magnitude, there 
is no ultimate solution to this problem, but it can be suffi- 
ciently mitigated to provide an overall high quality aerosol 
product. Moreover, the definition of AOT outliers is clearly 
application-dependent: for example, AOT enhancement near 
cloud edges may be a manifestation of 3-D radiative effects 
(Wen et al., 2006) or of complex aerosol-cloud interactions 
in the “twilight” zone (Koren et al., 2007), while being an 
obvious outlier for the climatology or air quality analysis. 

This paper described an improved cloud/snow screening 
technique in algorithm MAIAC which is implemented as part 
of the aerosol retrieval based on analysis of spectral residuals 
and spatial variability. Comparisons with AERONET aerosol 
measurements and a large-scale MODIS data analysis show 
strong suppression of aerosol optical depth outliers due to 
unresolved clouds and snow. At the same time, the devel- 
oped filter does not reduce the aerosol retrieval capability 
at high 1 km resolution in strongly inhomogeneous environ- 
ments, such as near centers of active fires. Despite signifi- 
cant improvement, the optical depth outliers in high spatial 
resolution data are and will remain a problem which needs 
to be addressed by the application-dependent specialized 
screening techniques. 
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